Graphene magnetoresistance in a parallel magnetic field: Spin polarization effect 
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We develop a theory for graphene magnetotransport in the presence of carrier spin polarization as 
induced, for example, by the application of an in-plane magnetic field (B) parallel to the 2D graphene 
layer. We predict a negative magnetoresistance a oc B 2 for intrinsic graphene, but for extrinsic 
graphene we find a non-monotonic magnetoresistance which is positive at lower magnetic fields 
(below the full spin-polarization) and negative at very high fields (above the full spin-polarization) . 
The conductivity of the minority spin band (— ) electrons does not vanish as the minority carrier 
density (n_) goes to zero. The residual conductivity of (— ) electrons at n_ =0 is unique to 
graphene. We discuss experimental implications of our theory. 
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Carrier transport properties of gated 2D graphene 
monolayers have been of active current interest [l], 0, 0, 
IH 0, @| , both for fundamental and technological rea- 
sons. From a fundamental perspective, understanding 
the mechanisms controlling graphene conductivity in ex- 
perimental samples is of obvious importance. Such an un- 
derstanding could lead to improvement in graphene mo- 
bilities, which in turn would help its eventual technologi- 
cal applications. The graphene chiral linear electron-hole 
carrier dispersion makes it particularly important that 
we understand the operative resistive scattering mech- 
anisms since the well-understood scattering theories of 
parabolic 2D semiconductor structures do not directly 
apply to graphene because of its peculiar band disper- 
sion described by the relativistic Dirac-Weyl equation for 
massless fermions. In this Letter we predict that the ap- 
plication of an external magnetic field parallel to the 2D 
graphene layer could lead to a novel magnetoresistance 
behavior in the low-density and high-mobility graphene 
samples, providing significant qualitative insight into its 
resistive scattering mechanisms. 

The physics underlying our predicted novel paral- 
lel field graphene magnetoresistance phenomenon is el- 
egantly simple: The applied field spin-polarizes the 
graphene carriers through the Zeeman effect, dramati- 
cally affecting the carrier density of states due to the lift- 
ing of the spin degeneracy, which then leads to a strong 
modification of the effective disorder scattering, thereby 
tuning the system resistivity. The key point here is that 
the parallel magnetic field couples only to the spin degree 
of freedom since, the 2D graphene layer being essentially 
a zero-thickness system, there is no magneto-orbital cou- 
pling to take into account as long as the applied magnetic 
field is strictly parallel to the 2D plane. This is different 
from the situation in the corresponding parabolic band 
2D semiconductor structures [9( where both Zeeman cou- 
pling and magneto-orbital coupling [Tol ] play significant 
roles in determining the parallel field magnetoresistance. 

Before describing the details of our theory and results, 



we start by providing a simple physical picture to moti- 
vate our predicted parallel-field-induced graphene mag- 
netoresistance, which could be either positive or neg- 
ative depending on the situation. Assuming that the 
only effect of the parallel magnetic field is to spin po- 
larize graphene carriers, we ask how such a continuously 
tunable (by changing B) carrier spin-polarization would 
affect the strength of disorder scattering by modifying 
the fundamental graphene parameters, namely, the Fermi 
wave vector (hp), and the Thomas-Fermi (TF) screening 
wave vector (qtf), both of which depend on the spin de- 
generacy factor g s . Since kp oc g s x ^ 2 and qtf oc gV 2 , 
we note that the Fermi (TF screening) wave vector in- 
creases (decreases) as the system becomes spin-polarized 
(with g s decreasing from 2 to 1) due to the application 
of the parallel field. Taking into account also the modifi- 
cation of the density of states (Dp) due to spin polariza- 
tion, we discuss below how the tuning of the carrier spin- 
polarization in gated graphene will affect the transport 
relaxation time, r, for various types of background 

disorder. 

Let us now consider the various types of resistive 
scattering potential V(q), which may be operational in 
graphene layers. For unscreened short-range white-noise 
disorder, V(q) is a constant independent of the wave vec- 
tor q, and as such, the changing spin-polarization would 
not affect the scattering potential, which gives rise to 
r oc kp . In general, an unscreened long-range disor- 
der potential would have V(q) ~ q~ n (with n = 1 be- 
ing the Coulomb disorder case for the charged impurity 
potential). Using q ~ kp, we conclude that the un- 
screened long-range disorder would effectively become 
weaker in the presence of a finite parallel field since 
spin-polarization would in general suppress the disorder 
strength by increasing kp, that is, we have r oc 
Finally, the 2D screened Coulomb disorder would have 
V(q) ~ (kp + qxF)^ 1 with kp (<?tf) increasing (decreas- 
ing) with the increasing spin polarization under the ap- 
plied parallel field, and as such screened Coulomb disor- 
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der could either increase or decrease in the presence of 
spin-polarization. Remembering that graphene conduc- 
tivity in the Boltzmann theory is given by a cx DpT, 
where r and Dp cx ^/g^ are respectively the scatter- 
ing time and the density of states at the Fermi level we 
conclude that if r cx kp 1 (i.e. if the disorder is inde- 
pendent of spin-polarization), then a cx g s , and it will 
always decrease (i.e. a positive magnetoresistance) in 
the presence of spin polarization. When the disorder is 
spin-polarization dependent (either through the modifi- 
cation of kp and/or (7tf), the magnetoresistance could 
be either positive or negative. In particular, for the un- 
screened long-range disorder potential, V(q) ~ we 
have a cx g\~ n - Note that for the bare Coulomb disorder 
(n = I) a remains a constant and does not depend on 
the spin polarization, implying no magnetoresistance. 

Below we develop a detailed quantitative theory for the 
graphene parallel-held magnetoresistance in the presence 
of disorder arising from random charged impurity centers. 
This is because the screened charged impurity disorder 
is almost universally [H, 0, 0, H, 0, @| considered to be 
the dominant scattering mechanism in limiting graphene 
carrier mobility in currently available samples. 

We hnd that the conductivity of intrinsic (i.e. undoped 
with Ep(B = 0) precisely at the Dirac point) graphene 
increases (i.e. negative magnetoresistance) with the ap- 
plied magnetic field. For extrinsic graphene we find a 
positive magnetoresistance for B < B s due to the mag- 
netic field induced changes of the screening properties 
of graphene, where B s is the field needed for full spin- 
polarization, and for B > B s we find a negative magne- 
toresistance due to the depopulation of electrons in the 
valence band. This negative magnetoresistance in extrin- 
sic graphene for B > B s , arising out of magnetic depop- 
ulation of the valence band, is the same qualitative effect 
which controls the negative magnetoresistance we find in 
intrinsic graphene. 

We use the Boltzmann transport theory including only 
the effect of resistive scattering by random charged im- 
purities 0| . The density of the random charged impu- 
rity centers is the only unknown parameter in our model, 
which sets the scale of the overall resistivity without af- 
fecting the p(B) dependence of interest in the problem. 
The Boltzmann conductivity of graphene at T = is 
given by a = e 2 v F Dpr(Ep)/2 where vp is the Fermi 
velocity, Dp = g s g v Ep /j 2 (7 = Hvp) is the density of 
states at Fermi energy (g s and g v are spin and valley 
degeneracy, respectively), r(Ep) = r is the transport re- 
laxation time, and Ep — jkp the Fermi energy. The 
relaxation time is calculated in the Boltzmann theory as 
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ing the background dielectric constant) is the graphene 
fine-structure constant, and e(q) = 1 + v(q)H(q) is the 
RPA dielectric function of the system, where v(q) = 
2ne 2 /Kq is the electron-electron 2D Coulomb interaction 
and Tl(q) the 2D irreducible finite wave vector polariz- 
abihty function. When the system is spin unpolarizcd 
with a 2D carrier density, n, the polarizabihty function 
becomes n(g) = DpP(q/kp) (HI, where 
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Then the conductivity becomes 

e 2 n 1 1 
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where I(x) is calculated to be 
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with 



sech \x) if x < 1 



When the parallel magnetic field is applied, the car- 
rier densities n± for spin up/down are not equal. Note 
that the total density n = n + + n_ is fixed by the ex- 
ternal gate. The spin-polarized densities themselves arc 
obtained from the relative shifts (i.e. the spin-splitting) 
in the spin up and down bands introduced by the Zeeman 
splitting associated with the applied field B. In this case 
the polarizabihty function becomes 



U(q) = D+P(q/k+) + DpP(q/k} 



(0) 



where Dp and kp are the density of states at Fermi en- 
ergy and the Fermi wave vector of spin up (down) state, 
respectively. Note that g v (—2 in the B — graphene 
case) is not lifted, but the spin degeneracy, by defini- 
tion, is lifted by the in-plane field B. The usual unpo- 
larized B = paramagnetic state has kp = kp = kp; 
n + = n_ = n/2. 

By defining the spin-polarization parameter given by 
£ = (n + -m_)/nwe can calculate the irreducible polariz- 
abihty function in the case of unequal population of n±. 
(Eqs. (l)-(3) above apply to the £ = unpolarized sit- 
uation.) In Fig. Q]we show the irreducible polarizabihty 
of the system for different spin-polarization parameter 
values. For unpolarized graphene (i.e. £ = 0) we have 
n± = n/2 and kp — ^2ixn± — ^fwn — kp. In this case 
the polarizabihty H(q) is constant up to q = 2kp and 
then increases with q. For fully spin polarized system (i.e. 
£ = 1) we have n + = n and n_ = 0; kp = \/2t\ 
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where fc F is the Fermi wave vector, m is the background and k F = 0. In this case H{q)/D F increases linearly up 
random charged impurity density and r s = e 2 /nj (n be- to q — 2k F with a slope (w/8)(kp/kp) — ttv2/16 and 
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FIG. 1: Total polarizability of the spin polarized system. 
Here £ = (n+ — n_)/n, Dp = gvEp/^j 2 , and ftF = \/2-Kn/g v 
is the Fermi wave vector of unpolarized system. 



has a kink at q — 2k p. For the partially polarized sys- 
tem, < £ < 1, H(q) has two distinct kink points at 
q = 2k$. 

Now we calculate the conductivity in the presence 
of a parallel magnetic held for screened Coulomb dis- 
order. In the presence of the magnetic field the to- 
tal conductivity can be expressed as a sum of conduc- 
tivities of spin up/down carriers, i.e., a = a + + <r_, 
where a± is the conductivity of the (±) spin subband, 
since the two spin channels are simply parallel trans- 
port channels as Coulomb scattering conserves the spin 
quantum number and therefore inter-spin-channel scat- 
tering is not allowed. The conductivities a± are given 
by <J ± = e 2 VpDpr(Ep)/2. In calculating the scattering 
times r(Ep) we have to use the two component polariz- 
ability given in Eq. ([6]). The total carrier resistivity p is 
defined by p = 1/a. 

First, we consider intrinsic (i.e. undoped or ungated, 
with n and Ep both being zero) graphene where the 
conduction band is empty and the valence band fully 
occupied at B = 0. Applying a parallel magnetic field 
to intrinsic graphene would cause a shift of the chem- 
ical potential from the Dirac point, and therefore the 
spin polarized (+) electrons induced from the valence 
band will start to occupy the empty conduction band, 
leaving holes created in the valence band. (The situ- 
ation is not dissimilar from what happens at T ^ 0.) 
The created intrinsic electron density, n, is the same as 
the intrinsic hole density, p, in the valence band, and 
the total density is proportional to the applied mag- 
netic field (B), i.e. n = p = (g v /16ir)(A/j) 2 , where 
A = g* psB (with g* is electron spin g-factor and pb the 
Bohr magneton) being the Zeeman energy. Since both 
carriers contribute to the conductivity and n — p, we 
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can calculate the total conductivity as a — ' 



h rii rj I(2r s ) ' 

Since n cx B 2 , the calculated conductivity of intrinsic 
graphene increases quadratically with the applied mag- 



netic field. Note that for this negative magnetoresistance 
phenomenon in intrinsic graphene as induced by the ap- 
plied parallel held, spin-polarization itself does not play 
any role - the physics is dominated by the field-induced 
creation of electron-hole pairs which can then carry cur- 
rent. 

In the case of extrinsic graphene (where Ep,n ^ 0) an 
applied magnetic field produces different numbers of spin 
up (+) and down (— ) electrons in the conduction band 
(or, for holes in the valence band) and we find 
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where p is the chemical potential, which is determined 
by conserving the total number of electrons, n = n + + 
n_. Then we have p{B) = jy/irn^l — S 2 /4, where 
6 = A/Ep and Ep — jkp — j^/im is the Fermi en- 
ergy at B = (unpolarized system). When B = 
(A = 0) we have p = Ep. We can relate the spin- 
polarization parameter £ to the applied magnetic field 
as £ = (5^/1 — <5 2 /4. At S = \/2 the system is fully spin 
polarized (i.e. £ = 1). Thus, the so-called saturation (or 
the spin-polarization) field for complete spin-polarization 
is given by B s = \^2Ep / (g* pb), which is smaller by a 
factor of \/2 than B s = 2Ep/(g* ps) for a 2D parabolic 
band dispersion. The conductivity of the partially polar- 
ized system can be calculated as 
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In Fig. [5] we show the calculated resistivity, p = 1/er, 
as a function of spin polarization for r s = 0.85 and 2.2 
which correspond to graphene on Si02 substrate and vac- 
uum. Note that we can represent the Fermi wave vectors 
as a function of the spin polarization, kp — kp^/l ± £. 
The total resistivity increases as the spin polarization (or 
magnetic field) increases because of overall suppression of 
screening. At £ = 1 (or B = B s ) we have 



I = 1 u r ° 

+ (l + 7rr s /4) 2 V l + 7rr s /4 
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and as £ — > 1 
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Thus, as £ — » 1, n_ — ► and the scattering time of 
spin (— ) electrons diverges as l/n_ because r± cx 1/J±. 
This means that the conductivity of spin (— ) electrons, 
cr_ cx ri—/I—, is finite at £ = 1 even though there is no 
spin down carrier. That is, we have a^(B s ) = 
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FIG. 2: (a) Calculated resistivity p(£)/p(0) as a function of 
the spin polarization for r s — 0.85 (solid lines) and 2.2 (dashed 
lines). Here p(0) is the resistance at B — 0. Inset shows the 
conductivities of each spin state. Note that the conductivity 
<t_ does not vanish at £ = 1 (or B = B 3 ). (b) Magnetore- 
sistance at B s , p(B s )/p(0), as a function of interaction pa- 
rameter r s . (Note that r s > 2.2 is not a physical regime.) 



which is independent of the interaction parameter r s . 
This is a consequence of vanishing density of states of 
graphene as n — * and the peculiarity of graphene 
screening properties, i.e. the contribution of interband 
transition to screening. We also show as an inset in Fig. 
[2] the calculated conductivities of each spin subband. As 
expected cr_ does not vanish as £ — > 1. The result in Fig. 
[2Ua) does not depend on the carrier density at B = 0. 
However, it is weakly dependent on the interaction pa- 
rameter r s . The maximum magnetoresistance at B a can 
be calculated as 



p(0) a(B s ) 
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where Jo = -f(2r s ), and is shown in Fig. [UJb). As 
r s — ► we have Iq,I+ — * t/4, and the ratio becomes 



P(B S )/P(P) 



1. As 



oo 7 — * 7r/64rf and i+ 



(16/7r 2 r 2 )/(4/7r), and the ratio becomes p(B s )/p(0) rs 
1.22. Thus we always have positive magnetoresistance in 
extrinsic graphene as the applied in-plane magnetic field 
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FIG. 3: Calculated resistivity p(B)/p(0) of graphene as a 
function of the in-plane magnetic field for r s — 0.85 and 2.2. 



increases from to B s . Note also that p(B s )/p(0) has a 
shallow maximum at r s rs 0.5. 

As the magnetic field increases beyond B s , i.e. B > 
B Sl the number density of spin (+) state is given by 



P 



n 1 
2^2 



8 1 
2 



1 



(13) 



where p is the induced hole density created in the valence 
band due to the magnetic field. In Fig. [3] we show the 
calculated resistivity of graphene as a function of mag- 
netic field. A very interesting feature is the sharp de- 
crease of resistivity for B > B Sl which is unexpected, 
when compared with ordinary 2D systems where the re- 
sistivity saturates (or increases ) for B > B s 0, E3- In 
graphene, however, for B ^> B s the calculated conduc- 
tivity increases as B 2 . This increasing conductivity for 
B > B s in extrinsic graphene has exactly the same phys- 
ical origin as the negative magnetoresistance we find for 
intrinsic graphene. 

Given the nontrivial (and surprising) magnetoresis- 
tive behavior of graphene we predict in a parallel mag- 
netic field, namely, negative magnetoresistance (<r(B) ~ 
B 2 ) in intrinsic graphene and positive magnetoresistance 
(p(B <C B s ) ~ B 2 ) as well as negative magnetoresis- 
tance (<r(B 3> B s ) ~ B 2 ) for extrinsic graphene, a 
question naturally arises about the experimental observ- 
ability of our predicted phenomenon. We believe that 
our predictions should be observable at low tempera- 
tures (< IK) in high-mobility graphene samples provided 
reasonable values of spin-polarization (e.g. £ = 0.5 or 
above) are achieved. Using g* — 2 for graphene we find 
A ps 0.1 2B meV where B is measured in tesla. This leads 
to B S (T) w 140v^ where h is the graphene carrier den- 
sity measured in unit of 10 10 cm" 2 . Thus, for n = 10 8 



cm 



B s 



14T, which should enable a direct verifica- 



tion of our theory close to the charge neutrality point. 
We anticipate g* to be enhanced by electron-electron in- 
teraction effects, perhaps as much as by a factor of 2 for 
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r s = 2.2, and therefore the necessary magnetic field val- 
ues may be a factor 2 lower for free-standing graphene 
in vacuum than on a substrate. It seems that a veri- 
fication (or falsification) of our predicted magnetoresis- 
tance behavior should be possible in high-mobility sus- 
pended graphene samples at low carrier densities [l3j |. 
This would help our understanding of the resistive mech- 
anisms and the nature of disorder in graphene environ- 
ment. 
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